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EFFICT OF YAPOS VELOCITY OR FILM CONDE SBATIO’ 
by 
James Ti. Webber 


(Submitted to the Department of Mechenical Engiveering on January 17, 
1955 in partial fulfiliment cf the requirements for the degree cf 
Master of Sctence in Mechanical Engineering.) 


In a vertical slate or tube condenser, the liquid condensate layer 
(film) flows downward under the influence of gravity forces. Laminar 
flow generally exists on the upper pert of the surface, but turbuient 
flow will exist on the lower parts if the condensation rate is out- 
ficiently high. If the vapor velocities are appreciable, the vapor 
exerts a shear stress on the liquid-vapor interface, which nay affect 
the heat transfer coefficients. 


The case cf laminar flow, under the influence of this vapor sheer 
stress, has beex analysed by Nusselt (1). By assuming the “woiverwal 
velocity distribution” of Prandtl-Nikuradse, Seban (10) obtained ane- 
lytical predictions for the turbulent fila, but neglected the effect 
of vapor velocities. This thesis extends the previous work by cevelon-~ 
ing equations whieh include the vapor sheer stress effect on a turbuioy 
film. Two eriteria for transition are proposed. Application of these 
eriteria persit the present resulta to be combined with those of Cusselt 
te give analytical heat transfer coefficients far the case of Laniner 
flow on the upper portion of a plate and turbulent flow on the iower 
portion. 


The results of the analysis are presented in the form of dimension- 
less plots, Figures 7 to 12, of heat transfer coefficieats versus Liquid 
Reynolds Sumber at the bottom of the plate and also versue plate length. 
Dimensionless shear stresses and liquid Prandtl] [hebers are parameters 
in these plots. The analytical results are shown to agree closely 
with Colburn's correlation equation (9) in the renge MN, = 2 to 5. 
However, considerable deviation exists in the very high” or very low WV 
range. The Prendtl Nwsbers covered by this analyais extend from 0.01 
to 10. 


Pr 


fn experimental set-up for investigeting the hydrodynamics of a 
falling liquid film is described. Preliminary results indicete that 






the analytical equations are able to predict film thick..38 quit) ac~ 
euretely for low Reynolds Nusbers. However, at higher Reynolds luebers, 
sane deviation was observed. It is concluded thet the Prandtl-iburedse 
veloeity distribution is generally adequate for the present purpose 

of predicting heat transfer ccefficients, but thet further stuily is 
necessary to verify this point. 


A proeedure utilizing a hot wire anemometer is described, which 


appears to hold promise in determining transition points in the i:quid 
film. 


Thesis Supervisor: Warren M. Rohsenow 
Title; Associate Professor of Mechanical Engineering 
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PAST I 


SECTION A - INTRODUCTION 

The influence of a vapor shearing strees on condensing film heat 
transfer coefficients has been shown by several investigators (1, 2, 
3, +) to be quite significant. However, this variable is often neglected 
in the interest of a manageable theoreticel analysis with the result 
that experimentally determined values of heat transfer coefficients 
have often differed from predicted values by severel fold. The devis- 
tion has been much more noticeable in the case cf vertical tube conden- 
sers where the condensate flows in film form on tke inside of the tubes. 
This seens reasonable when one considers the fact thet high vapor vel- 
ocities are more easily atteinable in this case than in the case of 
horizontal tube condensers or vertical condensers with the condensate 
in the shell. 


Although higher film condensate coefficients can be chtained with 





high vapor velocities it must be remembered that for vertical tube con-~ 
densera, the price of this increase in heat transfer ability ia an added 
pressure drop in the condensing section. In the case of steam conden= 
cers it is also to be noted that the cocling water side coefficient is 
usually controlling and little can be gaineéd in overall heat transfer 
coefficients by decreasing the condensate resistance to heat flow. How- 


ever, many industrial applications involving the condensation of organic 
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vapors are concerned with condensate coefficients as high as, or high- 
er than the cooling weter coefficients. It is in cases such as these 
that reductions in size and weight may be possible through the judicious 
use of this vapor sheer effect. 

A condensate film flowing vertically under the influence of grav- 
ity and a vapor shear stress lends itself remarkably well to theoret- 
ieal attack. It is the purpose of this thesis, therefore, to attempt 
to further the previous analyses by deriving expressions for average 
heat transfer coefficients which include beth laminar end turbulent 
condensate layers under the influence of a vapor shear streas. To sin- 
plify the resulting expressions, and at the same time give results of 
universel applicability, the case of a flat plate is considered. The 
equations may also be epplied to tubes es long as the thickness cf the 
condensate film is small compared to the tube radius, and as long as 


the effects of a pressure gradient may be neglected. 
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SECTION B - LAMINAR FIM 


In the folloving analysie the cage of a flat vertical plate, Fig- 
ure 1, is considered. Condensation is of a pure, saturated vapor with 
no noneondensible gases present. Assumptions include constant vall 
tenperature and constent fluid properties. The specific heat, *, » Gen~ 
sity, o, end thermal conductivity, k, may be determined et the arith- 
metic average film temperature. The viscosity, u, however, should be 
determined et the special average film temperature recommended by Col- 
burn (5), and defined by 


t =t - S/he >t (1) 


f 
The vapor shee stress, Tow » 18 consicered as being directly downward, 
and being constant along the length of the plate. 

The laminer enalysis essentially follows thet of Musselt (1), where 
the fila is in streamline motion and is under the influence of a vapor 


shear stress, Cail es show in Figure 1. 


Condensate 
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FICURE 1 
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By sssuming steady flow and neglecting mowentum changes, a force beli- 
ance for the shaded element shown gives, 
ad 
( --& (o -o,) (x, -«)- T.=0 (2) 
Ey v a) 
This gives the linear distribution of shear stress shown in Figure @. 


Sinee T = UT at x = 0, the well shear stress is given by 


Te = Be (o = a) *, * (ie (3) 


x 
FICHE 2 


The effect of the momentum crossing the element feces within the 
liquid is quite emall. To include the momentum effect at the liquid- 
vapor interface the term which should be edded to Rquation (2) is (1/¢, Max, df 
This represents the z direction velocity of the vapor crossing the li- 


quid-vapor interface, which is very neerly equal to the Liquid surface 
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a 


velocity, "2% end not the velocity of the vapor flcw in the main 
6 
vapor stream. This momentum term also is usually quite amall. 


Fron Newton's definition of viscosity, we have for leminar flow, 


av 
T ~~; ti a | 
&, ox (+) 
oe? 
Substituting in Equation (2) for e gives, 
6 i 
av, = B (op - py x, ~ xix +  bwax (5) 


This expression may be integrated across the film with Ti constant ; 


giving the velocity profile aa 


Pn 3a aod 
i... on ee ra a dhs as . 
v, = (o - ) Ca - 2/2) + bw x (6) 


The mass rate of flow per wit width of the vertical plate, [ , 


may be found by 


[ m (¢) ve ax (7) 
o] 
Carrying out this integration gives 1) in terms of film thickness, ae) 
4 2 
x ge. T % 
= «& a Sag wee B 
[eB |elo-0,) $+ 2; | (8) 


Equation (8) may be differentiated with reepect to x, to give @ rela- 
tionship between the change in mass flow rate per unit width and the 


change in film thickness, 


ca | 
Pe fee -eyh a -— | a, (3) 





a 


A linear temperature gredient may be assumed to exist within the 


condensate filma without introducing serious error. Keference (6) zives 
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~7— 
a complete treatment of the magnitudes i.ivolved in this assumption. 


Consider the control volume shown in Pigure 35. 


Temp U sat 
| 
t., | 
O 
x ~~ 
FIGURE 3 


If the enthalpy of iP entering through the upper surface of the con- 
trol volume is teken at a mean film temperature, it can be seen thet 
the [| leaving through the lower surface of the control volume will 
also be at this seme film temperature and hence will not enter into 
an energy balance. It is also evident that the enthalpy change of 

the increment of condensate, J! , will be equal to the enthalpy dif- 
ference between saturated vapor and a liquid at the mean Pilm tempera- 
ture. This may be designated be g and has been evaluated by Colburn, 
Millar, and vestwater (7) as he = De, + 3/8 ¢., St. On the basis of 
the previous discugsion an energy balance for the control voluse shown 
in Figure 5 may be made by equating the heat transfer to the change 


in enthalpy of the incrementel increase in condensate, De il , 
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faz = hy ae (10) 


Since the film considered is in laminar motion, the heat trens- 
fer coefficient may be taken as the thermal conductivity of the finid 


divided by the film thickness, 


k 
7 are * (11) 


Eliminating 9/A between Equations (10) and (11) and substituting for 
1" its equivalent from Equation (9), we have a differential expres- 
sion containing only dz and dx. as variebles, 


k At ss > 
— = cE e(p - e)x, + 








2 
wx 
6 





ax. (12) 


Equation (12) may be integrated from the top of the plate where x, = 0 


when z= 0, t& @ point Xye%e 





* 3 
Kot s Xo 9) . 
:” he ie de. ee _ (15) 
tg 
By rearranging end substituting the Prandtl number, i, ,® wk, the 


resulting expression is obtained, 


teat cp (3) - 45 @ 
Ne hy, (9 (I-¥) ( "S 5-4) (=) ow 


The following dimensionless lengths may be defined as 


. 


*o 
. = Wg) W> (15) 


Y 
es tas S& (x) 
N Pe hh’ p° \- = 
; (1-77) - 
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* 
Further define a dimensionless vapor shear stress, Te. a8 


2 tes . Fo fn 
“  a(p-Pe) (F)” 

J (17) 
Substituting these definitions in Equation (14) gives a simple and 
convenient expression for the variation of film thickness with length 
of the plate, 2. 

= (x,")* + 4/3 be (x,")? (18) 
The usefulness of the dimersionless quentities defined by Equa- 
tions (15), (16), and (17) may be more fully realized by deriving ex- 
pressions for the film Reynold‘s Number and mean heat transfer coef- 
ficient as follows. 


Define a film Reynold’s Number, Nae? bY 


(As) De . Ne (19) 
Ci 
By making use of Equation (8) for ‘, » we may write 


te. = a as = E 9(r- -f) Ae. - 5 (20) 


Rearranging the above expression, we have, 

3 “4. 

= as Xo 2 ds 1 Xo 
ae = a Re 8 hae coe 
ie -4 (C I Kr? 

7, 3 g 

Bis 
Now by substituting the previously defined values of os and _- 


Ne _ = 4 ae. 2. att 
(1- a 





(21) 


(22) 
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The original Nusselt relation gave only point values of the heat 
transfer coefficient and was difficult to apply in an actual design 
study as it involved a function of the condensate flow rate. ‘This 
Gifficulty can be avoided, however, if a mean cr average heat trene- 


fer coefficient is defined as, 


| hp 
. = ae (25) 


By substituting Equation (8) for [" in Equation (23) we have 


= bs diy ~fa- e Oo Ce. x 
Introducing the Prandtl Nuxber and rearran 


"3 ' 
ea ae es ee an 
# (5) > # Cy zat (AS) PL coy 43 ell 
(25) 


+ x” 
Making use of the dimensionless terms a. » X, , and (e » this reduces 


feo aie 3" ¢ 


a Ze 
(26) 
Equations (18), (22), and (26) completely define the results for 
the laminar region. For given values of the dimensionleas shear stress, 
ios » either the mean heat transfer coefficient group, or the Nie 
ean be plotted as a function of the dimensionless length of the plate, 


€ * 9 
2. This involves assuming values of x, end then solving for z from 
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Equation (18), He fros Bquation (22), ani he (x) frou Equation 
(26). This has been done for various A, and forms the basis for the 
curves to the left of the transition points in Pigures 7 to le. ‘The 
use of these curves in @ design case will we discussed in a later sec- 
tion, where the results of « turbulent film ae vell as the leminar 
film have been derived, and are plotted in terms of these fene varie 
ables. 
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SECTION C_- TRANSITION 


For a Liquid flowing in film form down @ vertical surface it seens 
reasonable to postulate that there will be transition from laminar flov 
to turbulent flow at a critical Reynolds number, in much the same muan- 
ner that transition occurs in « pipe full of fluid. When we speak of ze 
eritical Reynolds number it is to be remenbered that this is not a det - 
inite, distinct point at which transition will occur in all cases. 
Rather it represents some sort of a mean velue where the turbulence is 
sufficiently developed to materially affect the velocity profile within 
the film and where turbulent velocity fluctuations will be evident to 
the extent of affecting such characteristics of the filn as the heat 
transfer coefficient. 

Before proceeding with a discussion of transition, it becomes ne- 
cessary to define a Reynolds number to give a common besis of compari- 
eon. In the case under consideration it has become accepted to define 
the Reynolds number in terms of the masse rate of flow per unit cross~ 
sectional area, (w/S), the average fluid viscosity, uy, and en equive- 
lent hyéraulic diameter commonly taken as (48/Y)}, where Y is the wit 
width of the plate down which the fluid is flowing. In terms of these 


quantities, the Reynolds nuaber may be defined as, 


No. = Mee De a (72 1) a 7 (27) 


Although there has been little experimental work with the specific 
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aim of determining the transition point of a freely falling vertical 
film, certain points may be inferred from heat transfer data. In ref- 
erences (12) and (15) on the basis of sanewhet limited cata it is ree- 
ommended that a Reynolds number of 1900 be taken as the transition point 


for the cuse of zero vapor shear stress. It will be show later that 





considerable doubt exists as to this value, however from the present 
extent of knowledge no more precise value can be predictec. Therefore, 
for the purposes of the analysis, this criterion for transition will 
be adopted. One of the objectives of the experimental work underteken 
as part of this thesis will be to determine more adequate information 
on transition. 

Although we have adopted a critical Reynolds number of 1500 for 
zero vapor shear stress, there remains the considerable problem of con- 
sidering how a shear stress may modify this value. Carpenter (4) has 
suggested that under the influence of a vapor shear stresa, transition 
will ccewr at a given value of wall shear stress. This may be more 
readily visualized by considering the layer of fluid next to the wall. 
In laminar flow the sole effect of the outer layers on this inner layer 
will be a shear stress exerted on it. In effect, ali the inner layer 
“sees” of the outer layer is a tranemitted shear stress. If we first 
consider the case of sero vapor flow, but increasing liquid flow, this 
inner layer next to the wall will trip over to turbulent motion when 
the Reynolds number attains the value 1500. However, all the inner 
layer "knows" of this condition is a certain value of transmitted shear 


stress. 
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ow considering the case of a finite vapor shear stress, it is 

reasonable to assume that the leyer next to the wali will change to 
turbulent flow whenever this "critical" shear stress is exerted on it, 
whether by outer layers of lisuid or by an outer layer of vapor. Tak- 
ing the wall shear stress, fic » ®& the most convenient to which to ap- 
ply this criterion of critical shear stress, we may find its value from 
the known case of transition at a Reynolds number of 1900 with zero 
vapor shear stress. From Equation (3) ; Tocan be expressed in terms 


of x for Tr = 0, 
b: 4 _& — 
Th (p ey) x 


Similarly [* , end hence Reynolds number, + , may also be expressed 


in terms of x from Equation (&), 
te 3 
(4D) eran * 1800 = 4/3 Be (2 = a Jo) m, 


Eliminéting Xo between these two equations we ure able to arrive at an 


expression for the critical wall shear strese, designated by ( Te der nat? 
j9s0 (YA) 1% 
(Tegan © Ae (F- 9 | (y= ae (28) 


For the more general case where ie FO, we apply our postuistion 
that transition will still oceur at this value of a critical wall shear 
stress derived above. From equation (3) we have an expression for Ce 
in the presence of 4 vapor shear stress. By introducing the volue of 
wall shear stress obtained in Equation (23) and rearranging, the follow- 
ing expression for the dimensionless film thickness at transition can 


be obtained, 
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Substituting this value in Equation (22) gives an equation for the 


* 
transition Reynolés numbers in terms of the only varieble,  - 


WP Se - 
(Yemen * 1800 = 26(1-p,/o)/? Te. + 0.667(2 ~ p/p) (Te)? 


(50) 
This function is plotted in Figure 4 from a value of foe 2 to Te = 9.2. 


Inspection of Equation (29) will reveal that whenever e > Saaacaneea 
(i-p./p)"’ 
‘g 


the dimensionless film thickness, i will be negative at trensition 
Soviously this 48 impossible physically and represents a Limit to the 
applicability of the previovsly diseussed criterion ror traowition. 
In the higher ranges of vapor shear stress, we must, therefore, develop 
& new criterion. 

For the higher values of ‘e it is quite possible thet the film 


thickness at transition approaches a Limiting value, perhaps that of 





the laminar sublayer thickness which exists in ordinary pipe flcw. 


Carpenter (5) suggests the following magnitude for this thickness: 


% x e Jo Te z 
 ezen * PV? = 6 (32) 


Ordinarily a magnitude of 5 expresses the leminar subiayer thicimess. 
Actually an appropriate value is somewhere betveen 5 and perheps LO. 
The value of ». Is arbitrarily chosen. 

If Equation (71) is solved for To and this expreseion substituted 


in Equation (3), the following cubic equation relating film thichness 
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amd vapor shear stress is obtaines, 


(50"4)” + (x,",)° Tw oe (52) 
~ aye 


a 
may be found frem this equation through 


5 al 
For a given value of T- , Xo 


ge trial and errer process. With this value of x ‘ and the given value 
of "ga » the Reynolds number may be found through the use of Equation 
(22). The reeulting plot is shown in Figure 4 for values of cam 9.2. 
It can be seen that the assumption of transition at x” = 6 ie almost 
equivalent to the assuaption of transition at a Reynolds number of about 
70 to 30, there being Little variation of Peynolds mumber with vapor 
shear stress. 

Thie appears on first observation to be a somewhat stertingly low 
Reynolds muamber at which to expect transition, however, there is a log- 
jeal basis for its selection. There is very Littie good deta svailable 
to esteblish conditions under which transition occurs when there 1s a 
significant vapor shear stress. A single point can be obtained from 
the enthanol data of Carpenter (3) by plotting hie experimental points 
and noting the region of departure from the Nusselt laminar prediction. 
This transition point et « i ws 110 and Reynolds number of 96.5 agrees 
fairly well with the transition Reynolds number of 71.8 predicted under 
the present analysis. Tn terms of x." the transition data have the 
value of 6.9 which is close to the 6 of Equation (41). 

Aithough the criteria for transition ar portrayed by Figure 4 
have not been adequately substantiated, the resuits are nevertheless 


quite reasonable and are probably net sericusly in error. This entire 
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phenomenon of traneition as affected by vapor shear stress myst be in~ 
vestigated by obtaining much more experimental data. In spite of lack 
of adequate substantiation it seems reasonable to rT the transition 
erfteria as shown in Figure 4 for use in the present oi laa 
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SECTION D - TUFBULENT ANALYSIS 

When condensate flow rates become large, it seems reasonable to 
predict that turbulent flow will result, with considerable effect. on 
the heat transfer coefficients. It was observed experimentally by 
Kirkbride (2) that, in the case of vapor condensing on the outside of 
vertical tubes, measured heat transfer coefficients differed consider- 
ably from those predicted by the Nusselt laminar analysis. In these 
investigations there was evidently no appreciable vapor shear stress. 

Colburn (&) reviewed the results of Kirkbride and, by assuming 
that the sole resistance to heat transfer lay in the laminar sublayer, 
was able to derive heat transfer equations for a turbulent condensate 
or The Prandti-Nikuradse velocity distribution was assumed, and 
by selecting the thickness of the liminar sublayer at about x = 10, 
a satisfactory correlation of the data was obtained. This equation, 


as presented by Colburn (8) is 


Si J acecapd (t 7 aa 


It is to be noted that the Prandtl number enters here, while it did 
not in the laminar analysis. This does not seem ususuel when one con- 
: iders the correlation equations for heat transfer in pipes with no 
phase chenge. However, it should also be noted that this correlation 
was applied principally to Kirkbride*s data, with only amall variation 


of Nor in the range of 2 to 5. 
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The above analysis wes extended by Carpenter and Colburn (3,9) 
to include the effect of a vapor shear stress. In reference (3) heat 
transfer coefficients were investigated in a 5/8 ineh copper tube 9 
feet long. Special effort was made to accurately measure pressure drop, 
from which the magnitude of the vapor shear stress could be obtained. 
As the result of the experimental work and theoretical investigations, 
the following correlation equation for point values of heat trensfer 


coefficients vas obtained: 


V2 V2 


ie = 6.045 (New) TT. 
V2, 
k / (34) 
By integrating over the length of the tube, an equation for mean heat 


transfer coefficients was also obtained 


y /: 
bn = 0.065 (Ne) Te 
k 7” (35) 


By the introduction of the dimensionless vapor shear stress defined 


in Equation (17), this equation may reedily be converted into the form 


y Y Ya. 
1, a nx 
ha (E) ~ 0.065 (Np) (Tt) 
k Ad 
(26) 
The question of properly evalueting the vapor shear stress becomes 


somewhat complex when it is realized that, as condensation is taking 


place along the length of the tube, the vapor velocity is continually 
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decreasing. For convenience in making calculations, Carpenter and Col- 
burn reccommend taking a mean vapor shear stress, evaiueted from the 


usual plots of friction factor for one coaaponent flow at a vapor mass 





velocity, 
2 ° i/2 
» * & * & 
Ge = 
> 
(57) 
where G 


L 4s the v@alue of G at the top of the tube, and Go is its value 
at the bottom. It is felt that this represents somewhat of an oversin- 
Lification and that the friction factor should at least be evaluated 
from the limited data on concurrent two phase flow availeble (11). In 
reference (11) it is shown that there is considerable variation of 
friction factor from the smooth pipe data, mainly due to the fect that 
the liquid layer presents a rough surface to the flowing vapor. 

One point especially worthy of comment in connection with the work 
of Carpenter and Colburn is their consideration of the momentwa flux 
across the Liquid-vapor interface. In their derivation they evaiuate 
this tera as the product of the increment of condensate, d , and the 
main stream vapor velocity. This in effect assumes cone dimensional 
flow of the vapor. In view of the fect that a vapor shear stress is 
also considered, it seems sanewhat contradictory to assume, when eval- 
vating the momentum term, that a zero velocity gradient exists at the 
interface, and, when considering the vapor shear stress, to assume 4 


finite velocity gradient. It seems more realistic to visualize the 
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layer of vapor immediately next to the interface as flowing with eseen- 
tially the same velocity as the interface itself. In this case it can 
be reedily seen that such momentum terms will be small. It is through 
this line of reasoning that this momentwm term hes been neglected in 
the present analysis. 

(wite recently Seban (10) has applied analogy caleulations of the 
Prandtl-Karman type to a turbulent condensate iagee, essuming the 
"“wniversal velocity distribution” of Prandtl-Nikuredse. The improve- 
ment over previous analyses lies in the fact that the resistance to 
heat trensfer of the complete condensate layer is considered, rather 
than just the resistance of the laminer sublayer. Seban*s analytical 


results extend over the considerable range of H.. of 0.01 to 5 and 


Pr 
agree well with the empirical results of Colburn (&) in the range 

Ny, = 2to5. In the derivation Seban neglected the possible effects 
of a vapor shear stress. It is the purpose of the present anelysis 

to extend Seban’s results to include this effect and to complete ithe 
Story by obtaining mean heat trensfer coefficients whieh include a 
ieminar region at the top of the plate, followed by a turbulent range 
over the lower regions. 

Turbulence hes no effect on the force balance for a control vol- 
ume in the condensate film if we think of the shear stress as an apper~ 
ent TU  rether then the Newtonian viscous shear stress. Therefore, 
the resulting force balance equation will be identical with that for 
the leminar case, and is given by Equations (2) end (3). The same lin- 


ear distribution cf shear etress within the film will exist. 
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TT = (0 - o,) (x - *) + ce (2) 
9) 
‘i = . (p we Pry) *, a 3 ae (3) 


Since Newton's law of viscows shear will not hoid in the turbulent 
region, it is impossible to obtain analytically an expreseicn for the 
velocity distribution within the condensate film. This analysis will 
be based on the assucption that the velocity distribution will te that 


the Prandtl-Nikuradse formulation as given in reference (12) . 


For O<x <5 vo ax 
5 <x < ® v’ = + 3.05 + 5.0 In x” 
30 < x” vy’ = 5.5 + 2.5 ln x’ 
where: 
yy Sees oe =. jae 
yh oo I. 


(35) 


The apparent shear stress, i » for turbulent flow my be exproased 
os the sum cf the viscous shear stress ami the eddy diffusivity times 


the velocity gradient. 


~ 6 me onl. av = 
An { - * & ex ax ( 3%} 
wnere Cm is the eddy diffusivity. 
Rearrangin a dvr 
arranging, Sige (y+ é..) vr 
Substituting for T trom Bquetior (2) 


tel 3 s = We 
#13 (7 -f) (to-x) + Te] = (¥+ 8) & 
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Considering the region in the film where %0 < x x 


eseumed that << E,,. Solving for Cm JJ 3 


95 Car 9(1-7Fp) (Xe =<) 
ahi yr yy 
y/ {.- 


dx 
(40) 


The term dv/dx may be evaluated for the region under consideration 


(3<x < x) from the assumed universal velocity distribution, 


vo «= 5.5 + 2.5 nx 


v= Pe [5~ + 2.5 in ( /52)/ 


dv/ax = 2.5 | rte =z 
a Xx (41) 


The control volume for the energy balance is shown in Figure 5. 
Also shown is the type of temperature distribution to be expected for 


turbulent ?low. 
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For the shaded control volume it can be seen timt there will be no eu- 
taalpy ecnange for tie portion of condensate entering taroug:a the upper 
surfece and leaving tarough tne lower surface. Further, it is apparent 
thet the enthalpy chenge of the increment of condensate wili be essen- 
tially trat from saturated vapor to saturated liquid. Therefore, ac 


energy oelance will be of the fora, 


. % a 
a az Me ay 


As the left side of the control surface is moved across the film 
toward the wall, the above statements will hoid true except for the 
thin layer ext to the well where the temperature distribution experi- 
ences « sharp change in slope. Gince this layer is thin, the effect of 
subcooling may be ueglected for the turbulent film and q/A will equal 
(a/A).; or the heat transfer per unit area at the wall. Introducing 


this into the energy balance, we have 


af 
(+L? "we (N2) 


The equation for heat transfer across a turbulent film may be 
written in the form 
(Wa) « (a + €,) dt 
Pr ~ (43) 
where @ = = = thermal diffusion coefficient 
lao zs eddy diffusion coefficient 
Through the application of the heat - momentum transfer analogy that 


E, = Ey, , and knowing the velocity gradients throughout the film 
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from the assumed velocity distribution, we are able to obtain an expres- 
sion for the temperature gradients. By sucstituting Gy for €m in 


Bquation (40) and the velocity gradient from Equation (41), we obtain 
2%. — _9lI-"7p) (Xo-*) 





ate Ce a b 
ee ory y, ij 
di Xx 


Substituting this expression in Equation (45), the desired expression 
involving dt/dx and x as the only variables is obtained. 
g Jo es 3 /- tp) OF -x) 
( */A) | P P as 


+ 
~op¥ | New 454 [2 dex 





Solving for dt end introducing, for convenience sake, the following 


notation, B= = Ex 2 ie + 9 (I~ Lely) x 


nae Se ee 
PGR 1, te Bx 5 9 (l-P/) x 
Ne. a5V 5 eee 


(44) 


The condensate film may be divided into the following regions: 
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FICURE 6 


Equation (44) contains only two variables, t and x, since we have show 
that q/A will be constant across the film. This equation may then be 

integrated between the Limits “se to x, and x to %- After con- 
siderable algebraic manipulation and simplification, tne re cit of this 


integration across the film may be expressed az, 


Biv ty) op [9 (1-47) x : al it aM-l+)le oe 
M 


(FAs) ix 3 2M a) = <i 
fe M=- ve 
as Np. Xo 
66 ° 
eM eae” 

(45) 
_ [ _ } 
where M = <i °° ° 8 = ty 
aa. } x 
g ( P- 7) Xe | Ye 
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Equation (45) gives an expression for the temperature drop aecrose the 
highly turbulent core of the condensate film. By utilizing the results 
of a previous analysis by Yon Karmen (14) for the temperature difference 
between the wall and the cuter edge of the buffer layer (x = 30), we 
can obtain an expression for the totel temperature difference fram the 
liguid-vapor interface to the wall. Von Karman’s results may be put 


in tie forn 


(tx- tu) ae) VAC 9(I- H(I-Fp) Xe _ h Pr 
7 i ie . 5S Np. + 5 trl lt © Ne) 


(46) 
Adding Equations (45) and (46) te obtain the total temperature drop 


across the filn, 


ede / Cp je. 5 Net F ili 6H) t 
‘a ) 


jo ™M 60 | _ / 

25 iy |2M-1eV It Nee et i+ ee 
~ ae 10. TA £2 4 fy 10? = 10 

vs een "Vl ie 





(47) 
Let the right hand side of Equation (47) equal F,, then this ex- 


pression may be written 


at Yeas [3(1- PA) Xe = 
VA) M 


(47a) 
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For x => 9) 


, es r + 
[" = 7 (-6h + iu, + 25x, In x) 
For 5 <n < 30 


Po = PP (22.55 - 8.05 x + 


+ as 
5 x ln x 
aX . ) 


Differentiating with respect to x” 


m gives, 
For x" 7? W 


i 

i = fP(5.5 + 2.5 In x.) 
For 5 4a3< x 
pr 


= oe ss q sd 
dy,t 7P 3005 + 5 In .. ) 
Fer convenience define Fy as follows, 
+ aA + 
For — 7 30 Py = 5.5 + 2.5 in x, 


= -~ 5 e 
¥y = 3:05 + 5 In x, 


Then Equation (50) may be written as 


+ _ 
For x, > 30 are PUA 


For $<a,*< ye) 


For 5 <x < 30 al. PPE, 


Equetion 42 gives the results of the energy balance in the form 


(-B-) 


Oo 


ee. 


is may be written as 


ry = k, 2 sf 
A/, "72 at Jo 


430= 


(50) 


(52) 


(508 ) 










=e _——- =» "7 dla aL ay ota 1 
Soe te idl ho 
we —a) = «od a too pp 


_ Serena 
. SA 
wap abet eas 
SO lll 
“putea e # 
bende. | te 
rae me <"p 
” zz om # * ea 
ian} 


Se ee kee se ee ee 





- 

















a= 


Bo ec 


From Equation (47a) we have an expression for (q/A) as a function 


of film thickness, vhich may be rearreangec in the ~% 


(I- pyp\” 
g _ at Po (ke ) [ena (2)" y 
‘ i 
(53) 
Since it has been shown that (¢/A) = (o/A)., thie term may be eliminated 


between Equations (52) and (53). Substituting for from Equation 


xe 
(50a), we obtain an expression containing, as the only variebles, ~~ 


+ 
ax, » and dz. 3 Ys 


ae lama es GE) Gel * 


Thies may be integrated between the limits x * to x. end z. to Z, where 


ot 


oe and Z, indicate the values at the transition from laminar to turbu- 
lent flow. The reeulting equation may be put into the convenient form 
xt Vp 
fi, Les = a ma? 


vs Gh |- 77/7) 
'/3 
ied 7) (54) 


Equation (54) now gives the desired expression for the distance down 





(Z-Z) At = 


the plate, 2, as a function of film thickness. In solving for this 

cag 
distance, the value of x, at the transition point must first be solved 
for from either Equation (29) or (32) and then the distance to the trans- 


ition point, z,, found from Equation (18). The distance from the 











~Le- 


Seen a me eet Ar Lea) TAN) Veh 
mica tt gitervomee 120) ban (8) 


: i 7 AOR eR ae ae 7 


xh ~) Fa = -&h 7 ie 
We Ae a} Seale 
ree Pt teak a at Sd dled ld detined , 
ret oma sav end? ll ba emaen it wit 
2. 7 a he 





oe 


LE eT jeer bevien car enemy ue (48) Gee 

1h oN nntnine at tance AED He ie © ae ‘ 
ois seal aoa A 
ee ee 


- * 


‘ 
5 





Boe 


transition point to the point under consideration may then be calculated 
from Zquetion (5+) by evelusting the integral. This must be done nunmer- 
ically. 

While point values of the heat transfer coefficient es given in 
Equation (43) are of academic interest, they do not give answers readily 
usable by the designer. To predict the heat transfer characteristics 
ef a condenser or to enable rational design, it becomes necessary to 


predict mean heat trensfer coefficients which can be applied to the 





ndenser 8 a whole rether than et just one point. Accordingly it 
is possible to define this mean heat transfer coefficient in the sane 


manner es was done for the wholly laminar filn. 


| 


h, = | = DR hes 


LY (tx, - Eu) ; L(ts- Cu) 


(55) 
It should be noted at this point that the substitution of q = w . 
ef Dee implicit in Equation (55) neglects the effect of subccooling 


of the condensate film. As has been shown previously this ie very near- 
ly true for the turbulent film, but it is not true for the leminar film. 
Since, in the present case, we are considering both a laminar and s 
turbulent film, this will introduce some error in the results. There 
will be no error in the leminar curves provided an equivalent Dee a6 
defined in Equation (10) is substituted for he, The error involved 

in neglecting subcooling may be somewhat evident in the region just 
following transition, tut will become less importent as more and more 


of the film is in turbulent motion. In most of the cases encountered 
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the leminar film and transition region wili extend only a few inches 
from the top of the plate und the neglecting of the subcooliing in this 
email area should not appreciably affect the overall results for plates 
whose lengths may approach several feet- It is also felt that the un- 
certainties of actually determining the transition point will introduce 
errors of at least as greet a magnitude as the subcooling error, and 
therefore this latter point should not be belabored. 

Since we have assumed a constant temperature difference, the de- 


nominator of Equation (55) may be broken into two parts. 


h,, = i he 


ee a | ee 
(tx. te) acer (ti, - ti} (L-2:) (56) 


The term (t = t.) z, applies to the laminer portion of the condensate 
fs) 
film, and the term (t ~ tL = z,) applies to the turbulent portion. 
O 


Fron Fquetions (16) end (18), 


, We ; 
Z(tx, - t) = (I-%7p) New hy, (F) [ (as) + 4 c* (424) | 


4 Cp 
(57) 
Equation (54) gives the expression for the turbulent length. Sub- 
stituting these two values into Equation (5%) ytelds the final equation 
for the mean heat transfer coefficient in terms of functions of the 


film thickness, 
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(58) 

The Reynolds Sunber at the bottom of the plate is of some interest 
as it is a direct measure of the amount of vapor that has been condensed. 
This Ree nm t and may be readily evaluated from Equations (49) 

For < VA 30 


- e ‘ ¢ * 
Noe w «256 + 12 Xo + 10 Kor la Ke. 


% =< 
For 5 < Xo L 
+ 4. a 


Hae = 50.2 - 32.2 , + 20 ... in Bo. (59) 


Bquations (54), (58), end (59) represent the resulta of the tur- 
bulent analysis in much the sase way that Equations (15), (26), and (22) 
represent the results of the laminar analysis. The seme forus of plot- 
ting are therefore used, the turbulent plots merely being continuations 
of the leminar curves. A description of the mmerical evaluations of 


these equations is contained in the Appendix A. 
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SECEION E - DISCUSSION OF ANALYTICAL RESULTS 


The results of the analysis presented in this thesis are contained 
in Figures 7 through le. These curves ere besed on a rational approach, 
end show the effect of a vapor shear stress on heat trensfer in laminer 


The 


& 


and turbulent film condensation for a wide range cf Be = end ie? 

use of these curves will be illustrated later by a sample problea. 
It is significant thet the plotted values of mean or average heat 

trensfer coefficient include the effect of a leminar condensate film 

on the upper part of the plate and a turbulent film on the lower regions, 

if one exists. To use these plots it is not necessary to first inves- 

tigate whether or not turbulent flow exists, this determination having 

already been made in the preparation of the plots. The values given 

are wiiversally applicable and are not dependent on the seonetry of 

the systea or the fluid properties except as they enter into the deter- 

mination of the dimensionless quantities needed to enter the curves. 
Inspection of Figures { through LY reveal several interesting points 

concerning the ve of film condensation. First, the significant 

effect of a vapor shear stress ie apparent. This effect is more pro- 

nounced in the lower Re e T88, giving en increeze in heat transfer 

coefficients of several folé in sane casea. At Higher Rae the effect is 

less. The effect on h,, of the exact loeation of the transition point, 

es influenced by the vapor shear stress, becomes less significant as 


the value of a. > at the tottom of the surface increases and is practically 
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insignificant for N. » greater than 10,000. 

The effect of turbulence in the film is illustrated by the change 
in slope following the transition points. It can be seen thet the tur- 
bulence effect on the heat trensfer coefficients is greater for increas- 
ing Ne, 2 and for increasing i, The Reynolds Number effect might be 
expected, as larger Ap ‘ will result in a greater degree of turbulence 
and hence a greater portion of neat transfer by the mechanism of tur- 
bulent mixing. The Prandtl Number effect is best seen by comparing 


the heat trensfer plots for the very low N, = 0.01, representative of 


Pr 
Liquid metals, and the plot for the high N, = 1C. The curves show that 


Pr 
turbulence hes a greater role in increasing the heat transfer character- 
istics for the higher qe* In fact the low Prendtl nuwsber curves for 
the turbulent portion actually lie below an extension of the laminar 
curves. This may be explained once it is realized that for these lig- 
ulds, the thermal conductivity is so large that it becomes the principal 
meens of heat transfer, rether then turbulent mixing. Since the essumed 
velocity distribution results in a thicker condensate film then would 
exist for leminar flow, it can be seen that this thicker film offers 
& greater resistence to heat flow through conduction than would the 
laminar film. 

It would be desirable to compere the present analysis directly 
with experimentel results. However, experimentel data on film conden- 
sation in which the vapor shear stress has been eccurately determined 
was not readily available. There does exist a correlation equation 


based on considerable experimental data, which has already been mentioned, 
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and is given by Equation (36). In order to campare the analytical 
results with this correlation equation, the curves in Figures 8, 10, 
end 12 heve been cross plotted at a Reyncids Mumber of 5600 and are 
presented as curves of mean heat transfer coefficient versus 1 ie . 
Theee plots are represented by the solid lines in Figure 15 for various 
Prandtl Nuwabers. The correlation equation in the form given by Iqua- 
tien (36) hae also been plotted, and is show by the detted lines. As 


can be seen, the agreement is quite good for HN, = 2 to 5, which lies 


Pr 
within the range of 2 to 5 of the liquids actually tested and on which 
the correlation equation is based. There is considerable deviation cut- 


side this renge, as evidenced by the curves for #. = 10 and 0.01. This 


Pr 
indicates that caution should be employed when using the correlation 
equation outside the range for which it was proposed. 

Figure 13 also graphically illustrates the considerable effect 
thet a vapor shear stress may have on the heat transfer coefficients. 

For exemple, the h. group for _ s 10 and (ie 7 50 is about 1.00, or 
four times the value of 0.25 given for T. = 0. 

In the force belance given by Fquation (2) the only pressure grad- 
ident in the « direction considered wes that due to the height of the 
vapor and is represented by the = fr tera. While this assuaption 
was veliad for the case of a flat plate in an infinite atmosphere which 
we considered, it may lose its validity when a vertical pipe is consid- 
ered. For pipe flow, the existance of a vapor shear stress is intimately 


connected with a pressure gradient in the z direction. To include this 
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=) 5m 
effect, it is neeessacy to rewrite the force balance, Equation (2) in 
the form | 
ae ae. ye = 
T ~ (7 - GE) (%o~X)- T. = 0 


(60) 


To evaluate ap » consider a control volume containing the vapor 


poe 
in the core of a vertical pipe, as shown in Figure 14. 








Condensate 
film 
FIGURE 14 
Considering the vapor flow as one dimensional and including the change 
in momentum, the pressure gradient may be evaluated from a force bal~ 


aace on the control volume as, 
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(61) 
where u_ represents the vapor velocity. 
The weight rate of flow of the vapor, VW, =P, Ux a which may 


be differentiated to give, 


ee _. for AW 
iz {z 


The weight rate of flow of the liquid condensate, w = 2 x ror which 


may also be differentiated to give, 
de pr all 


da dz 
Since ice = - 3 » We may solve between these two equations for the 
2 


velocity gradient, 
a ee Om 


ia fe ga i 
Substituting this expression in Equation (61) gives the pressure grad-~ 


tent as 


apa 2A § 2. Welt  _ at. 
daz Jo OF Ve ode re 
(62) 
Previously, in the case of the flat plate, we had taken the pres- 
sure gradient = (e/g) p,- However, for the ease of the vertical pipe 


condenser we may take sn equivalent vapor density 


(63) 
This equivalent vapor density, which takes into seceount the pressure 


gradient for a vertical pipe, contains two terms which are dependent 


a in 
on T,, the radius of the condensate layer interface. For a given tube, 
the equations presented in the enalysis can be evaluated in a manner 
analogous to that for the flat plate. Since at and u will vary down 
the length of the tube, these must be determined from the given conditions 
at each step in the stepwise numerical integration. It is apparent 
thet the solutions so obtained will apply to only the particuler gean- 
etry for which they were calculated and for the particular vapor vel- 
ocity which wes considered. In this procedure, the advantage of the 
universality which was obtained for the flat plate will be lost. 

In Fquation (63) the terms containing the vepor shear stress and 
the vapor velocity are subtractive end may quite possibly largely can- 
cel out. If the net megmitude cf Equation (63) is small compared to 
the condensate density, p, it is then possible to apply the results 
obtained for the flat plate. 

It is worthwhile to investigate the magnitude of the pressure 
gradient correction for e specific case. Consider a 1 inch diemeter 
tube, in which water is condensing. Saturated vapor temperature will 
be taken 2s 120° F., and the tube wail temperature as 60° F. To take 
an extreme case, assume (hs = 50. 

Meen film temperature, fron Samition (1), 

t, = 120 ~ 3/4 (60) = 75° F. 
Then 

= 1.000 X 107” rt”/see 

ge 1.46 x 107" ft. 


Since 4. << / 
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ieee a point in the tube where x = 50. Then x. = 500 end 
LA] on 
At 75° F., the Prandtl Nusber of water is 6.29. Proceeding in the man- 
ner outlined in the section on calculation methods, FP, may be evaluated. 
F, 2 57 6 
Then taking k of the condensate es .357 BTU/nr-ft °F, 
h = 2660 BrU/mr -rt~ - Pp. 
Substituting h, he = 1026 BTU/Ib, and at = 6O°F., 


Af = 2.59 Ib/nr - #t° 
ys 


Taking smooth pipe friction factors for illustrative purposes, 
f= .005 at a vepor Reynold's Number of 64,500. Since Cweet & 
we may solve for the vapor velocity, u, 

u = 1090 ft/sec 


(Tho by = fu UD 2a 64,500 


AZ” 
Substituting our computed values in Equation (63), we are able to 


a 
evaluate Py » 


+ 2.34 ~ 21.02 nl 19.48 1b/tt? 


wey * 
Py 205 
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Although p,* is sbout 1/3 of the condensate density, ¢, this is 
not the magnitude of the error involved by neglecting this term. Se- 
ferring to Equation (60), by continuing our present specific example, 
we can evaluate the maximum error in the shear streas, i? which will 
eecur when x = 0, or T =als - ince .~ = 50, .< 50 X 1.46 X 10°" 


rt. al -0073 Lt. 


9 
(A = i. Te (= 7? - dB ) Xo 
To « 0455 + (62.4 + 19.48) (.0073) 


(Tom +455 + -59T = 2.052 iy/ee® 
Ry neglecting the pressure gradient, 
G3 incor, © °H55 + «55 = «910 Ib/tt® 

The final error which is introduced in this specific case is of 
the order of 15.5%. These calculations are included, not to exhaust- 
ively investigate the effect of the presswre gradient term, but rather 
to give en order of magnitude of the error which might be expected. » 
The values chosen in the problem are felt to be representative except 
for the vapor velocity of 1090 ft/sec. This high value wee deliberately 
pieked to indicate en upper limit to the error involved. In any prac- 
tical steam condenser, the vapor velocity will undoubtedly be much lower, 
with lower vapor sheer stresses. It may be concluded that in many cases, 
the eifect of the pressure gradient may be neglected, end the enalyti- 
cal results for the flat plate may be applied tc vertical tube conden- 
sers with confidence. 


In all of the calculations for Figures 7 through 12, the quantity 
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(1 - p_/p) was taken as wnity; @.g., p_/p emell. It was also assumed 
that Tr and (t ~ t_) were wiiform elong the plate. The equations 
in the atibents itis be readily evaluated secomting for any value of 
(1 - e,/e) and for eny proserived variations of (t_ - t_) or of Tr 
along the length. The effects cf these variations _ functions of 
the particular geometry and are hard to generalize. 

Of the meny asguaptions underlying this analysis, most are sound 
and heve considerable experimental evidence to buck wo their validity. 
The two most subject to question are the assumed velocity distribution 
and the criteria for transition. It will be the purpose cof the second 
section of this thesis to investigate the hydrodynamics of a falling 
vertical Liquid film under the influence of an appreciable shear stress 
at the interface. These investigations exe preliminary only and will 
merely point out the direction to be followed in a more detailed exper- 


imental program. 
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SECTION F - SAMPLE DESIGN PROBLEMS 


Application of the results of the analysis is conveniently 
made with the use of the dimensionlees plots of Figures 7 through 
12. Since mean or overall heat transfer coefficients are dealt 
with throughout, no point by point or stepwise calculations ere 
necessary. The actual methods involved in a desisn celculation 
are illustrated ty the following problems. 

The two problems selected ere felt to be representative of 
those which might be met in industrial practice. The operating 
conditions and tube sizes chosen correspond to two actual cases 
where measurements have been taken. The calculetions should 
therefore serve the additional purpose of cunpering the accuracy 
of the predictions made by the present amalysis. The first prob- 
lem utilizes tne initial date for water rm 191 of Carpenter's 
experimental program, reference (3). The second is based on an 
actual industrial design as reported by wurster, reference (15). 

The problem is posed in the following mammer. For a given 
condensate weight rate of flow per tube, it is desired to deter- 
mine the required length of tube to condense the amount speci- 
fied, with no uncondensed vapor passing out of the bottom of the 
tube. In addition, it is required to compute the mean condensing 
side heat transfer coefficients. 


In both cases the condensing vapor is stews at essentially 
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atmospheric pressure. The wall temperature has been teken as 
the arithmetic average of the wall temperetures reported in 


eech of the references. 


Ho. 1 No. 2 
Tube Cherecteristics: 5/8" ,14 BYG 1", 1 7BYG 
Inside Diameter 0459 834 in. 
Cross-Sectional Area .00115 00426 re” 
Tube wall temperature, t, 162°F 182 .8° F 
Desired condensate flow rate, w 50.1 201 lb/hr 
Vapor mass velocity at top of tube,G, 45,600 (47,100 == Ib/nr-£t* 
Since G,°9, Eq. (47) gives Ge 25, 300 27, 400 1o/nr-re* 
Specific volume of saturated vapor at 14.7 pia, v=26.8 ft? /U 
Mean vepor velocity, Un 188 .2 204 ft/sec 
Condensate flow rate per unit perimeter 
at bottom of tube, |" =w/Y 418 870 lb/hr -ft 
Average condensate flow rate I'm 209 435 lv /hr =ft 
Since water is the fluid being condensed, the ratio of surface 
tensions, s/s = 1. 
The density of saturated liquid at 212° F., 9 = 59.7 lb/ft? 
Then = 3.5 7.28 et /ur 
The viscosity of steam at 212° F., y= -0303 lb/ft br 


Cm a 
The Feynolés Number of the vapor, a 29,700 66 ,G0O 


Utilizing the two conpon 
(11), f 


nt flow friction factor data in reference 
-0063 O72 





z 
ZU 
Vapor shear stress ie ft tn -1296 174 1b/et" 









att - ‘ 


<n.amay uae ne agua Davee ea eae J rmepients 
Scr eehertyginoet sees 
wait a 

a 
| oferta, 6) eee 

at ih -~. Be Smet eee 





~ 















hes. ed | oe nemoeny 
vySar ~ ar a i 
sal \ead (oe _~ .\ewr er | ’ 


ewido ee paul im tbat lap 
a a ee PS ICY oe os me 
(OWT 0 teee yalee VAP deeper Dee 


“ 
















=o ie ao -? ry 
ee a 
Ten ors Me. a ot et o 
ew ar uw e- a Sr wey eommmanes egenee 
abe. & Lider air ,omemmieom mhée hfest aie 0) votes 
ee | 
Vale Ta, Sp ie alee oe a 
wih . as 7 . = an _ 
nt ibd i - We pes 


i wae Th a oe a ee a 


wren oem — © ere 


ae ‘ 
“omnes ve ou. “ah 2.7 ee Gee eer 








~53- 


iO.) Ose 
Mean Film Temperature, from Eq. (1), te 17%.5° 190° 
a i «% 0 
Kinematic viseosity et t., 9X 10° .%2xK 10° ft/sec 
} 1 
(e)” 190 X10"? .74l x 107" ft 
jig = nA oe 27.4 49.5 
A * Pi 
9(/-7+) (ey 
ak am 2 
Viscosity of condensate at t., 1! 154 % 107 .680 X 107? 1b-see/tt* 
Cerdensate Reynclds Nuzver at 
hottom of tube, Nn. = + 1915. LAZO 
Pranétl Number of Condensste at t,, N,. 2.26 2.02 
Pater Fig. 10 for W,, = 1 with N, 
% Pi > 3 me 
and Ss , 5 get hn (5) 4] 47 
K \3 
Enter Fig. le for Re 2 10 with Roe 
a *yH/3 
and Te 3 get fe (E) 7% 33 
) } hn (. Is ‘7 
Interpolate to ectual Ki, , ne ( ; -} AG 38 
Therwueal Conductivity of condensate Bot, 
at t., k - 386 WO pn iet OF 
Mean heat transfer coefficient, h_ 2250 2470 
% ayl/3 t 
Z = ‘he he (+) L160 3350 
Since z = 4¥zat Cp a} la 
Nepr hs V )-f4/ 7) 
solve for z 2 ( 5.20 10.4 ft 
Actuel 2, as determined by experiment, 
(references (3) and (15) 5.6 11.7 ft 


By comparing caleulated and actual tube lengths, it can be seen 
thet the analytiesi results will predict the performance cf a ecndenser 
quite closely. Probably the greatest source of error lies in the selec- 


tion of the correct average vapor shear 
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stress end tube wall tempereture. [If it was desired to predict 
performance within very close tolerances, the numerical inte- 
gration procedure described in Appendix A could be reworked with 
the known variations of 4 +t and te included in a stepwise approx- 
imation. However, this would involve a rather tedious and lengthy 
set cf calculations end, for most eases, it is felt that the con- 
stant average values as used above will provide results of suf - 


ficient aceuracy. 
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PART II 


SECTION A - LITRERATULE SURVEY 


Part I of this thesis is concerned with the developaent of ana- 
lytical expressions predicting the heat trensfer coefficients of a 
falling film condenser. Use nas been made of the heat momentum trans- 
fer enalogy, so that the validity of the predicted heat trensfer coef- 
ficients is dependent upon a knowledge of the hydrodynamics of a falling 
film. Assumptions have been made regarding the behavior of such a fila, 
based on the best knowledge available today. It is the purpose cf 
thie section, therefore, to give a brief review of the literature cov- 


ering the subject of falling liquid films. 





unread and Badger (13) discuss the onset of turbulence in a felling 

liquid film. Their conelusion, based on the work of other investigators, 
is that the critical Reynolds Number is in the range of 1600. ‘They 
also cite date of Jakob and Erk (19) as showing that transition mey 
cecur at lower B > in the presence cf high values of vapor velocity. 

In reference (17), Hebbard and Badger measured film coefficients 
of heat transfer on a 1 inch 6.D. vertical tube condenser le feet long. 
Reynolds Nusbers cf the Pilm covered the range 500 to 3000. Results 
are consistently higher than those predicted by the Kusselt Laminar 
equations, and generally lie in ea line parallel to the predicted line. 
io measurements of the vapor shear stress were made, and it is feit 
by this author that the discrepencies may quite possibly be explained 


by this effeet. No information as to the onset of turbulent flow can 
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be deduced from these results. 

Kirkbride, reference (2), presents heat transfer data of several 
investigators on Liquid filme flowing down vertical walls. This data 
covers the Reynolds Number range of 700 to 20,000. By the marked dif- 
ference in slope cf the correlating line from the Russelt laminar pre- 
dietion, turbulent Plow of the film ia definitely indicated. However, 
there is little data for the lower nN. e” 2° it, ig not possible to direct-~- 
ly infer the Location of the transition point. The best correlating 
line through the data intersects the theoretical Nusselt prediction 
at a Noe of about 1100. The vapor shear stress was not considered. 

In addition to considering heat tranefer coefficients, Kirkbride 
(2) presents data on measured film thicknesses for MN. lees then 2000. 
Excellent agreement with predicted thicknesses is obtained for the lower 
Ni 


Re 
ere given as the explanation for the discrepancies. 


» however, deviations are observed for the higher Bae? surface waves 


Cooper, Drew, and McAdams (15) present isothermal data of eix 
observers dealing with Liquide flowing down the wetted smooth surface 
of vertical towers or over flat plates at slight inclinations. Film 
thicknesses were messured either by a micrometer, visually cbheerving 
contact with the film, or by suddenly stopping the flow and measuring 
the quantity that drains out of the test section. In a plot of the 
Fanning friction factor vereus Noe? scatter cf the date from the lemi- 
nar equation starts at a i, 2 of about 1000, with the best correlating 
curve indicating transition at abeut 2100. Data were not available 


for determining the possible effects of gas velocity on the Liquid flow. 


~Vc= 


eer samt? cowl sented og 
Lemire, & atu) legend amb Sean ,(f) eeestos . stryahue 

wed ofl vali laniree ome ques «LID bogul @ eeeesgidere 
-Shh ee ee DO OR te eer eet GAlsoeee GBP send 
75 ek FL Oe aA wh! eke Per ITT Mar SD meDTe as OmoeeT 
we er antkes B49 TRA Of B01 ad Wo tt feta eerbid 
Si ee ee 
SRO Pees ear Lamtog euetiRy er he muteoel er SeweE et 

ee Fie Aer eet ermetornt bred «at dpdenur eats 
See ee ee ke se 
<RoSHUT tee Thee eee Peel uolwshteneny Sf schtnhh Am 
cena t etunena 
oem eam. SAE WE Ker eeOD Oe ROLITEIIND eererae Ie 
‘ i eer Aa nt mieretem att a owt) be 
nod Ve ated comet eomeevy (AI) eenlhed hee sree) yee 
























mkrt seO0O MEAG) TOME Te AetAy Gat Swe oc viwel tectbew Ye 
WER Chm OO EYE MOE Semen Vrer mY 
i ee ee ee. Y 
FD OT eth Teme Kemet ee Fe Rey Mee Sey CPE romney A 
SO Ek A Ee Lee ST eT SNE 
IAC Ha eT SS ym) Teele ee ee Nee NEES ie 
bt Nee | i mm Ui rF aap To arte efvines, ae) guftetmemet wt 











~58= 

Considerable experimental work has been done on non-condensing 
heat exchangers employing vertical tubes, in which the fluid flows in 
film form down the tube wells, references (20) to (26). Bays in ref- 
erence (Zé) has plotted heat transfer coefficients versus H_° A break 
in the curve at a MN of 1500 or 1400 may possibly indicate transition 
of the film from leminar to turbulent flow, however, Bays recommends 
a value of 1000. Leavitt's data (25) for water extends from aN, of 
1000 up, with no indication of transition as low as 1000, although 
considerable scatter of the data is observed below Noe « 1700. The 
date. of Fentress (24) covers the Reynolds Number range of 850 to 20,000, 
with no break observed which might indicate transition. In ail of theese 
references, no vapor shear stress was present and no attempts were made 
to measure film thicknesses. 

Dukler in reference (27) has used a film thickness measuring tech- 
nique employing a capacitance measurement. This device determines film 
thickness by the difference in dielectric constants of air and water 
between two plates of a capacitor. One of the advantages of the tech-~- 
nique is that it everages out surface waves and irregularities. All 
measurements were made on the film flowing down a vertical, flat piate. 
Film thickness readings agree well with the Nusselt prediction in the 
lower HN. range and with the prediction from the "universal velocity 
distribution” in the higher = range. Asp the technique is somewhat 
difficult to apply in a round tube, no investigations of a vapor shear 


stress effect were undertaken. For the flat plate with no vapor shear 
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stress, trensition is indicated at a N.. of 1000. 

It is apparent that considerable variation exists in the results 
of the many investigators cited here. On the tasis of the works cov- 
ered in this survey, no difinite conclusions can be reached as to the 
onget of turbulence in a Liquid film flowing downward under the influ- 
enee of grevity. Yurthermmore, no investigations have been made to de- 
termine the effects of an appreciable vapor shear stress. More knowl- 
edge in this direction is urgently needed. This, therefore, will be 
the direction to be taken in the experimental part of the thesis. 
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The purpowe of the experimental progran is to investigate the 
hydrodynamics of a liquid film flowing downward elong a vertical sur- 
face wider the influence of a shear strese exerted om the cuter face 
of the liquid. In order to cbtain quantitative results it is necese 
sacy to be able to directly measure the following guantities: 

(1) The liquid flow rate 

(2) Tne Liquid temperature 

(3) The magnitude cf the interfacial shear stresa 

(4) The liquid fila thickness 

(5) Im addition, it would be desiraplie to be able to obtain 

measurementa which would give quentitative information as 
te friction factors for two component flow 

Although a two dimensional set-up, with the liquid flowing down 
a flat plate, was considered, the necessity of accurately determining 
the interfacial shear stress led to the choice of a vertical tube for 
the test section. By introducing the liquid along the tube wall, it 
fiews down the valis in film form. A flow of air in the core of the 
tube sinmlates the vapor flow in @ vertical tuve condenser and supplies 
the shear streas on the outer face of the liquid. By asesuning no 
pressure drop across the liquid film, this sheer stresa can be easily 
measured by pressure drop readings along the length of the tube. 

The test section is shown in Figure 15 and ia a one inch dismsneter 


tube ayproximetely § feet long. Pressure taps are installed at one 
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foot interveis along ite length. In order to prevent ore ing over of 
the liquid at high flow rates, & sacoth convergent nozzle is placed 
mt the top of the tube. Surrounding the nozzle is a settling tank 
to damp out eny turbulence or fluctuations which might exist in the 
Liquid ag it is discharged from the supply system. The nozzle ects 
as a weir in the settling tank to aveniy distribute the liquid along 
the walle of the tube. The entire test section, nozzle, and settling 
tank are constructed of lucite to permit visual gteervation of flow 
cherecteristics. a 

Six probe ports are installed in the test section. These are 
arrenged in two groups of three each, spaced evenly around the cir- 
ecumference of the tube. Details of the probe perts are shown in 
Figure 16. Film thickness readings are made with the modified depth 
gage micrometer, showm in place. Metal probe plugs » Whese inner ends 
are flush with the test section wall, are installed to provide elec- 
trical continuity from the micrometer probe through the liquid film. 


The micrometer probe is e standard depth gage micrometer with a 





pointed extension added bo ‘the shaft. The shaft and the extension, 
except for the extreme tip, are insulated. This probe is shown in 
Figure 17. To obtein film thiekness readings, the probe is inserted 
in tha port, with the shoulders of the micrometer hela firmly on the 
blocks provided. A sero reading is obtained by advancing the micro- 
meter until it touches the fax wall of the test section. The probe 
is retracted end a sensitive Chometer commected between the end of 


the micrazeter and one of the adjacent port plugs. The micrometer 


~S0— 


wero griews doeresy of vets al wtp of) geola ADeyoetal foot 
heenlg ol aloo: Jepgneveds dooce « ,aetwo welt yin 7a biapil or 
dnd gollivee «2 bt eieeme ai) gelkeeewert «pee ade Yo Get at? Je 
el) of fetee gio Goldy sookimetee( tO serelered we tuo peed oo 
atee eizeoe 26% .atrype Chega wey met) fequmuieeto et 22 ae Dips 
pools blietl ese afodiarecd Clover o) dae pedis a oO ee ee 
pristine See alone Jhvows Seer ee a oe eur to ali oe 
oul? ke mmiterseds tamaty Time w etm te herertees eee Set 
9 Pt eee 
ow saad? .culsoes tens oo) «ch Setieeeel we atu eortg €De 
-1i8 aod Meee elesve beeeg> ine erwie Se ager ont at beg 
ob wads ote ahcg oe al ie oot i eee 
Udgee A)Uihow mil Sie eld ome Spinone eemmiY WEE «OS waryee 
ahte see omy lg Mey Leow eo of ede Jodeeoreis eee 
-99)9 milvewg of Setiooe fee Liew eosees feed eds Ay cele oe 
ALN) DLopLi wet ORO? along AERESER eA Mey GleKIeOU dames 
8 Ale Talweroly Mem Cee Pests 6 6) ory wOMmIoals eT 

~edihomtee sdb cos Mote ot? | ml of Dwkde uloeEe hetasag 
nh won wi sdocg eit? tla Ow (ghey eter anf ol cee 
Dacenaed aS Gees SE) ,xpediows eomamecee wll) wiatdes of «TL ome 
aay ce gia D Bled mrenotois eat Yo wedded: ec dilv . frog al GS 
-125ke <r goltervie qi dewtadds «t exter, ovee 4 .kebivew amole 
og «ff .apevows fees wit te Liew ot add eadquet fl itz: wow 
& bee at eweitet pefoeces weeC ryliieoee & fae beforevot af 
eerste a oat fay Seek be to moe Bee cetera sale 

















a 


























-63- 





ee a Et. 





zd IY, Vi 
GN 


Cai 
LU, 











VW 


( 
aN 






Sw 
th +4 (> 
t} a) {le QO, t q 
ob ee 
i fogs 
=} 36 
2, 
re 
rf) 
a? 
1¥ : 
La 2 
C 
oe C4 


fly 
@) 
re) 
< 
at 








(o> 


eqoId JO VworO THY 


 — ae 


= 
a 


ee ee 


= 
= - 
S 








alii 


is then advanced until electrical continuity is established by the 
probe tip touching the outer edge of the liquid film. This eontin- 
uity is easily detected by the ohmeeter. The insulation on the shat't 
cf the probe prevents completion of the cireuit by the film on the 
neer wall of the tubc, through which the probe extends. 

Pressure readings in the test secticm were made with a standard 
U tube manometer vith air in the tubing between the manometer end the 
pressure tap. Surface tension forces at the interface between the 
liquid on the walls of the test section and the air in the tubing, 
pilus the tendancy of the Liquid te be draw from the test section in- 
to the tubing when changing the pressure level, made these readings 
somevhat difficult. Uncertainties of measured pressure levels vere 
only evident at low air flow rates, when the surface tension forces 
were of a comparable magnitude to the pressure forces. Modifications 
have been made to correct this by using a single menameter tube, with 
the connecting tubing completely filled with liquid. This requires 
that o zero reading be taken with no flow in the test section, but 
eliminates all interfaces in the pressure measuring apparatus. How- 
ever, all data presented here were made with the original syeter. 
Pressure readings plotted es straight lines for all but the very Low 
air flow rates, and appear to be fairly reliable. 

The test section discharges to a closed separating tank, whose 
purpose is to separate the liquid from the air. The air suction line 


is connected to the top of the tank through a standard pipe Tee inside 
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the tack. The air flow must traverse several sharp cie:.ges i.. ddreetion, 
thereby separatiig a:y liquid droplets widen might be entrained. A 
Liquid recirculating pump is connected to tne bottan of the seperating 
tank through a flexible cose. 

_ So attempt was made to measure tre sir flow rete in tuese prelim- 
inary investigations. Suen measuremecits could be easily made by instal- 
ling a rotameter or some positive displaceme:t meesuring cevice in tase 
air suction line. It is felt that trese measurements, which would en- 
able deterninatio: of frictien factors, would be nignhly desirable for 
future experime::tal work. 

The Liquid flow rate is measured by an orifice installed in tse 
Supply line. The pressure Grop across this orifice is measured by & 
water over mercury mancseter. This system was calibrated in place by 
time - veig:.t measurements. Calibration data and curves are given in 
Appendix B. Liguid temperature readirgs are obtained ol a thermometer 
placed in the settling tenk at the top of the test section. 

Auxiliary equipment consists of « pump to recireulate the liquid 
through the system, and necessary pipixg. For the preliminery experi- 
mental progrem included in this thesis, vater wae the only liquid tested 
end therefore, connections were mee to a water main es a sowree. A 
line sketch of the equipment is shown in Figure 18. Figure 19 is a 
photograph of the set-up, showing the test section assembly and sepa- 


rating tank on the right, with the micrometer probe in piacc. 
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Figure 19 


Experimental Apparatus 
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Since it wes anticipated that vibrations would quite possibly be 
critical in effecting the steadiness of the flow, flexible hoses were 
instelled in the piping from the separating tank ta the pump and in 
the supply piping just before the settiing tank. The test section was 


thus isolated from vibrations in the piping system and pump. 
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SECTION C - FILM THICKNESS 


Pilm thickness measurements have been made in the manner described 
in Section B, Part TI. Since we are dealing with @ round tube rather 
than a flat plate, it would be worthwhile to investigate the equations 
applicable to this case and compare them with those derived for the 
flat plate in earlier sections. 

Consider the control volume, whose cross section is shown shaded 


in Figure 20. This control volume is an annular ring eround the tube. 


Vapor 


Condensate 


PIGURE ©2O 
Since there is no condensation, momentwe terms may be neglected, 


anc the foree balance for the control volume is given as 
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Tatrd2 + {2 dan (y~re) E ; fit(r-K)dz -T,. aK dz =0 





(c+) 
After some simplification and rearrangewent, this reduces to 
3 2 iat @ is _tdp (=e 
ala r + To} a @ r (65) 


Since monentun terms do not enter, it is possible to use the ex- 


pression for db/Je of Equation (61) with the second term = 0 


a. i Oh Ee 
dp. a / 


ri 
- (66) 
Substituting this in Equation (U5), we get 
9 (r=) ee 
- 95 (/-%,) T- F Arr: 
(67) 


In solving for dp/dz from conditions in the vapor core of the tube 
and then epplying this dp/déz to the equations for shear stress in the 
liquid film, the assumption has been made that no change in the pres- 
Suré gradient oecurs in a ra@fiel direction. ‘hie assumption is quite 
good, and almost exact fcr thin liquid files. 

In order to compare the distribution of shear stress given by 
Equation (67) with thet previously derived for the flat plate, it is 
necessary to substitute ry = ~ 3 xX, end rin roo Xe These substitu- 
tions are identically true from the geometry and merely serve to shift 


the zero axis from the tube centerline to the tube Wall. Making the 
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substitutions and simplifying, we have 


Yo(%o—*) (Xo ~ Xx) ex 
T= (7 1) wer ih | eT, ee Me 
(65) 

If the liquid film is thin, ro > x and Lio? Xp For this econdi- 
tien considerable simplification may be made, resulting in the foilow- 
ing shear stress distribution. 

= (7-7) (xx) + Te 
(69) 
This expression is identical to Equation (2), which was derived for 
the flat plate, end forms the basis for the development of the predicted 
heat transfer coefficients. 

In the experimentel progrem, the meximum film thickness measured 
was of the order of .025 inches. This represents 5% of the tube radius. 
By substituting these vaiues in Equation (68), it is easily shown that 
the maximaun error which is introduced by the simplifying assumption 
thet the film is thin is only 5.%. This maxinmuwa error will eccur when 
considering the wall shear stress and will only be apparent when the 


vapor shear stress forces are large compared to the gravity forces. 





In this case, the film thicknesses ere much smaller than .0c5 inches. 
Having verified that the flat plate equations for shear streas 
are epproximately true for the case of a vertical tube, it is now pos- 
sible to utilize the flat plate equations to predict the film thick- 
ness as a fimction of i - Por leminar flow, Equation (22) applies. 


# 
For a given Reynolds Number, by asswiing various values of X and 
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solving for age we are able to plot _ as; @ function of %.. 
Similarly, for the case of turbulent flow in the liquid film, for a 
given Reynolds Nunber, « trial and error solution of Equation (49) 
will yield a value of x, . From the definition of x Beuation (33), 
the wall shear stress can be found in terms of the film thickness, xo° 
By equating this expression for vell shear stress to that fran the 
foree balance, Equation (3), a formula for 7 in terms of To 46 
obtained. This has been done for both laminar and turbulent films at 
tne Reynolds Nunbers for which film thickness readings were obtained, 
and ferms the basis for the theoretical curves shown on the plots of 
measured data. 

The procedure used in cbteining measured values of film thickness 
is as foliows. Liquid was recirculatec with the pwap, end the flow 
rate adjusted to give the Feynolds Number desired. The pwap wes used 
on all of these runs, rather than water directly from the water main 
es it was found that fluctuations in city water pressure made it impos- 
sible to obtain steady flow rates. With the wate: Plow rate constant, 
the throttling vaive on the air suction line wea aijusted to give var- 
ious levels of pressure drop through the test section. At each of 
these values of pressure drop, film thickness readings were taken with 
the micrometer probe, following the procedure described in Section B, 
Part II. Before each series of runs, film thickness readings were 
taken through ail three of the probe ports in the group being vorked 


with to insure that film thicimesses were equal around the circumference 
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of the tube and that the flow was evenly distributed. The raw deta 
and computations for these runs ere included in Appendix C and are 
shown plotted in Figures 21 through 27. 

The static pressure level at points in the test section was plot- 
ted egeainst distance, and a straight line curve drawn through these 
points. The slope of this curve was then teken as the pressure grad- 
ient in the z-direction. This pressure gradient cowld then be easily 


* 
converted to the vapor shear stress,  - , trom the formula 


a * ®t ® tb 
‘(ae * * 
(70) 


where _ is the radius of the tube. 
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USSTON OF FIM THICKNESS MEASUREMENTS 





Inspection of tae plotted results of the film thickness measurc- 


ments, Figures 21 to &7, reveals several interesting fects. Firat, 





the shapes of the curves of x " versus i’ correspond ¢losely with 
the shapes of the thecretical curves. In general, the values of meas- 
ured film thickness lie beiow the predicted values, this deviation be- 
coning more noticeable at higher liguid Reynolds Numbers. The measured 
values for e = 0 are consistently higher then predicted. At high 
values of a ,» experimental values of x” agree well with theoreti- 
eal values. There appears to be rather sharp change in curvature of 
the measured values at values of iA in the neighborhood of 5 to 10d. 

For + O, visual observations revealed that rather lerge sur- 
face weves were present on the interface of the liquid layer. These 
vaves made film thickness measurements, using the micrometer probe and 
electrical continuity technique, somewhat uncertain. The ohmneeter 
registered alternute breaking end making of ccntaect between these sur- 
face vaves and the probe. In taking the readings, an attempt was made 
to odtein a mean wave height reading, but this procedure is adnittedly 
@ crude one, as surface tension forees tended to maintain contect be- 
tween the probe tip and the liquid surface. Therefore, it is not sur- 
prising that these readings were high, it veing believed that thick- 
nesses, as measured, were actually somewhat nearer the crests of the 
waves than the mean height. 


” 
For all values of E greater than about 5, the large, siow moving 
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surface waves were replaced by such smaller and more rapidly moving 
irregularities. Indications on the chnmeter geve quite sharp cut-off 
points for contact of the probe with the surface. It is felt that film 
idescds daktanes Pee We Hider Wie Oe 1. are cents 46 Tithin 
-0005 inches. This is based on the sharp cut-off points and the con- 
sistency of the readings on several rus under identical conditions. 

The fact that measured values of ." lie below the predicted cur- 
ves leads to the conelusion that the velocity distribution is actually 
something other than that assumed, Although the theoretical turbulent 
curves lie above the leminer curves, it is believed that, for the higher 
vaines of e the film is essentially in turbulent motion. This can- 
clusion is based on the appearence of the film and on some prelininary 
results obtained with a hot wire anemameter in the liquid film. The 
fact that the velocity distribution may be something other than the 
universal distribution of Prandtl-Nikuradse is not altogether surpris- 
ing. The wee of this velocity distribution is certainly an over exten- 
Sion of their work, which wes derived from measurements of full pipe 
flow. Nowever, the actusl values of film thickness are eat least with- 
in a few percent of the predicted results. 

It is believed that the change in curvature of the measured values 
of film thickness is due to the demping out of the large surface waves. 


This corresponds to visual observations. The general cheracter of the 





flow, as observed through the transparent test section was as follows: 
* 
For Tt. = O, the flow was meinily cle@r with relatively lerge and reg- 


ular surface waves. The frequency of these waves was low and fairly 
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constant, being of the order of 5 cycles per second. As the air flow 
was increased, ‘these large surface waves were demped owt until, at 
around ths 5, the flow appeared quite clear and transparent. For 
yet higher values of t » very suall surface waves were formed, with 
@ relatively high frequency. It is estimated that the height of these 
waves wag less ‘then .0005 inches. It is felt that these amall weves 
are directly connected with the onset of turbulence in the liquid film. 


It is to be noted in this respect, thet these small waves could be 





eny air flow at all, but required very high liquid 
Reynolds Nyabers. No measurements of this point were obtainec, the 
flow rate being beyond the capacity of the liguid flow menaneter. It 
is estimated that liquid Reynolds Numbers of the order of magnitude 
of 6,000 to 8,000 are required to attain this condition. 

Gn the beeis of the film thickness meesurenents, no definite in- 
formation could be deduced as to the onset of turbulence, as effected 
by @ vapor shear stress. However, the equations predicting film thick- 
ness proved to be approximately correct. It may be concluded fran this 
information that the assumed velocity distribution is adequate for the 
purpose of predicting heat transfer coefficients. It may also be con- 
cluded that the equations derived for a flat plate will be equally 
applicable te vertical tubes, as long as the film thickness is small 


compared to the tube radius. 
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SECTION D - PRELIMIZARY REQULES OF HOT WIRE ANEMOMETER 


In Part I te criteria for deveraining the transition point cf 
the liquid Pilm from laminar to turbulent flow were proposed. These 
criteria, although based on same experimental work, are far from being 
verified, especialiy in the region of epprecieable vapor shear stresses. 
In facet, no experimental work hes been found thet can shed light on 
tris phenomenon of transition, eas affected by a vaper shear stress. 
This section will diseuss the preliminary work for determining the tran- 
Bition point by means of e@ hot wire anemometer. No quantitative infor- 
mation has been obtained, but work along these lines is curremtly in 
progress and should become available in the near future, reference (16). 
The instrument is a standard hot wire anemometer, manufactured 
vy the Flow Corporation cf Cambridge, Measachusetits. It is being used 
in ccumressor blede cascade studies in the Sas Turbine Laboratory at. 
M. I. T. The element is a .GO03 inch diameter tungsten wire, supported 
ty a long shaft. By inserting the shaft threugh one of the probe porta 
en the test section, the wire can be introduced into the liquid film 
on the fax side of the test section. With currents of approximetely 
150 milliempa, good indications can be obtained om an osellicecope. With 
the relatively high current, the wire burns cut almost instantly if ex- 
posed to air. This nes been the souree of seme difficulty, and it may 
be necensary to modify the wire element or lower the current to pre- 
vent this. 


At Low liquid Reynolds Mmubers, the flow appears to be quite 
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susceptible to unsteadiness ceused by vibrations. Preliminery runs 
were made while e large air compressor wes rummiing in the seme build- 
ing, approximately 100 feet avay. Very regular velocity fluctuations 
were observed at low Keynolds Musbers, which died out as the liquid 
flow was increased. ‘The runs were then repeated at night, when the 
compressor was not running, and these Tliuctuations were absent. 

For the case of no shear stress, Laminar flow in the liquid film 
existed up to Reynolds umber of 6,000 to 5,000. A trip wire was in- 
serted around the nozzle throat at the entrance to the test section. 
With this trip wire, cecssional turbulent fluctuations were observed 
et the mexismm flow rete that the liquid flow manometer could measure. 
Yt ia believed thet stable laminar flow exists up to Reynolds Numbers 
of 6,000. This observetion has at present not been verified. 

With air flowing in the test section, « consistently turbulent 
flow cowld be attained with lower Liquid Reynolds Numbers. However, 
as wes evident fram the film thickness measurements, the film is ex- 
tremely thin for the regions of appreciable vapor shear stress, and 
difficulty wes encountered in keeping the hot wire completely in the 
water file. It way be necessary to use a more viscous fluid in order 
to get film thicknesses of sufficient magnitude with which to work. 

On the basis of the very sketchy work done with this technique, 
Little of a quentitative nature can be coneluded. However, the tech- 
nique appears to be a powerful one. Turbulent indications arc quite 
warked on the oscilloscope trace, although a lower beckground noise 


Level woulé be desireable to eneble more exact determination of the 
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transition frem completely laminar conditions. It is felt that, once 
the minor probleas have been eolved, the hot wire anemameter will yield 


the answer to the transition question. 
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SECTION & - RECOMMENDATIONS FOR FUTURE EXPORDGINTAL Wonk 


In view of the present lack of knowledge concerning the behevior 
of e falling liquid film from a fluid dynamics point of view, it is 
recamsended that the first experimental efforts be directed along this 
chennel. 

Specifically, it is felt that it would be desirable to fully in- 
vestigate the trensition criteria. This should be done with et least 
three liquids of widely varying viscosity. Consideration should be 
given to varying surface tension also. This reeomendation is based 
on the fact that, in this type of fluid flow, these forces may quite 
possibiy be the only significant omnes in the radial direction. It is 
worthy of coement that a surface tension ratio was the correlating 
factor in the friction factors presented in reference (11). 

In working with other fluids, it would be desireble to verify 
the film thickness measurements presented in this thesis. By gather-~ 
ing sufficient information, equations more adequately describing the 
velocity Gistribution in the turbulent regions should be developed. 

Much more information eoneering friction fectorsa in two phase 
eocurrent flow is needed. As previously pointed out, this can be 
eesily obtained fron the present set-up by adding a meter to measure 
air Tlow rates. 

After the hydrodynasics of the liquid film have been sufficiently 
defined, the analysis should be modified to include any new or better 


information thet may neve been discovered. Finally the analysis should 
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be directly checked by investigation of en experimental condenser, in 
which provisions should be meade to determine point values snd mean 


values of the condensing side heat trensfer coefficients. 
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APPEADIX A_- NUMERICAL EVALUATION OF ANALYTIC ESMATION 
Laminer Portion 


For this portion of the condensate film, the significant Equations 


are (15), (22), and (26). 


22M +4 Oey (18) 
wine § 1s) «gt OT 
(\-*/¢) - (22) 
be (Z)?. & (ty ee os) 
tm. (F) - $ +4 Ved - 


For given values of T eesume various values of x." lower than the 
transition point, which is given by either Bquation (29) or (32). For 
the re Moe ° x » both 2 and oe Tp) may be evaluated 
directly. 7 ( 5 } may now be found from the given values and 
the computed values of z . 
Turbulent Portion 

For computation purposes, Equation (53) may be rearranged into 


convenient form 


hn ( iH Nee 


ok NS [m- 22) +20 





the more 


kK 
(55a) 
The evaluation of this equation entails e numerical integration pro- 


cedure enc is best followed by listing the steps involved. 
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* 
Step 1 - Fer given vaiues of (2 ad Noe? aseume various values 
% 
of x greater than that at the transition point. 
Step 2 - Evaluate M from Equetion (45a). 


+ og * 
Step 3 - Evaluate =X,  &, + Te l= o/o) . at this 


point in the canputations it was necessary to asewre that 
(1 - p/p) = 1. 

Step 4 - Calculate F,, from Equation (47). 

Step 5 - Calculate F, from Equation -- 

Step 6 ~ Form the product, FF, ka ‘ 

Step 7 - Plot this product as a function of x,» 

Step 3 - At various values of x", take the area under this curve 
with @ planimeter or by other suitable graphical meens. 
Thie area is taken from the transition point, Xa to 


the selected valves of x x! 


"V5 
Step 9 - Caleulate (2-2) = + mr pi | A d Xe 
Np (I- yp Xe ’ 

Step 10 ~- Add to the result of step 9 the distance to the trans 
ition point, ane found in the leminar calculations. 
Step 11 ~ Caleulate Noe from Equation (59). 


) 
Step le - How be (1) r Nre 
hk \Y pal 


kere 2. is the result of step 10. 
hn (2)" : 
Step 15 - Piot e fe) versus £ end A, ° There will result a 
* 
double series of curves vith T. and i as peraneters. 
The results of the above calculations are contained in Table I and are 


Shown plotted in Figures 7 through 12. 
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TABLE I 


ANALYTIC VALUES OF BEAT TRANSFER COEFF IC Tawrs 





= OL 
= 0, e. - 14 ,900 
= : 
a. Ka 
30 2k 580 
100 50,500 
200 265 , 700 
300 5 30, 500 
400 857,700 
606 1,662,900 
800 2,626,500 
1000 5,668 ,500 
i. = 2.5, z. = 7420 
~~ oe 20,560 
100 72 , 460 
200 248 860 
280 502,620 
4OO 817,820 
600 1,595,420 
800 2 , 528 , 620 
1000 5, 57%, 620 
Tei = 5, a © 215 
— | 
100 68,015 
200 237,615 
5300 461,615 
400 786 , 415 
600 1,539,215 
800 2,442,015 
1000 3,463,615 


2, We 

5 ote 
12, 744 
20,464 
28 , 54nb 
46 ,2h4 
62,844 
80 , 344 


2, 302 

5, 54h 
12, Ti 
20,464 
28 , 54s 
46,244 
62 , 844 
Bo, S44 


2, 302 

5 Sa 
12, 74 
20, 464 
26,544 
46, 2h4 
62,844 
80 , Bly 


09565 
206887 
204796 
-035357 
.055e8 
202761 
02 993 
202204 


01125 

-O7651 
05121 
4071 
2035490 
02899 
202435 
262 


01250 

O3151 
05 363 
04249 
03650 
-03004 
025735 
602 534 
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T, = 10 e. = 80.9 e)/s 
$e . 
s. - Re Bae (5) 
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40 730,882 28 , 544 95905 
690 1,445 , 281 4G 2h 03200 
500 2, 508,081 62,644 02725 
L000 5,276,081 60, 8444 02468 
vO = 20 7 «= 61.5 
30 4,862 1,48 “23514 
50 13,662 2,2 L635 
190 52,062 5 She 1065 
200 190,862 12, 74 -06677 
300 396 ,862 20, 454 «05156 
400 856,062 26 , 544 oO43551 
600 1,301,662 46 ahh 203553 
800 2,083,662 62,844 035016 
LOCO 2,972,462 60,844 02720 
a) eT 
50 3,240 1,125 Te 
50 9, 640 2, 302 22388 
100 33,040 5, 54s 1457 
200 Leh O80 12,744 03348 
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APPEUDIX B - CALIBRATION GF WATER FLOW ORIFICE 


Bec. 
76.6 
63.0 
71.6 
T9 4 
160.5 
171.0 
Lis .5 
GO.4 
245 .3 
79.4 
69.7 
58 69 
52.6 
54.5 
Lg 8 


Lb/sec 


0653 
01589 
61398 
L260 
00623 
0585 


TABLE If 
Teup F GEM 

DB 471 
13 1.142 
ro 1.007 
73 905 
15 ig 
13 A421 
19 2640 
19 197 
19 12gh 
79 2909 
79 1.054 
19 1.222 
19 1.370 
79 1.329 
79 a) 


inches Hg 


2.08 
12.05 
GAB 
7.63 
1.95 
1.65 
3678 
5 86 

77 
7 48 
9 98 
13.97 
17655 
16.47 
19.65 
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APPENDIX C - BATA AND CALCULATIONS POR FILM VRICKNESS 


Theoretical Predictions of Film Thickness 
laminar - The laminar predictions result directly from Equation 


(22), which may be rearranged in the form 


4 *\F 
as. _Nre = 3 (4) 
yy 

*% 
For a given Rae » the plot is made by assuning values of X and solving 
push 
for lw . 
Turbulent - The turbulent predictions result fron the following 


procedure: First, consider Equation (49), which relates Nand ,, 


For = > BW 





Wy a a_./* = (-64 + 3 x, t 2.5 x ln x, ) 


+ 
Fora’ ¢x%, < 0 


ak ia Mpe/* = (12.55 - 6.05 x 4 5 x 


“an x" 
7? o) 


oO 


For a given Nae . 


Now, by eliminating ‘To between the jefiuition of _ Sa tquation (38), 


this equetion is solved by trial and error for x 


we 
and the force balance, Equation (3), Ow may be expressed in terms 


1 (Og) «. x 
tes (1- *-/) Ue : 


é x.* wha. 
@) an 
Xo Xo ’ 


With the value of x," determined by the Ne » and by acesuming various 


*% * 
values of x, T. @an be solved for. 
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The above procedures for both the leminar and turbulent films will 
yield the predieted film thicknesses shown plotted in Figures 21 to 27. 
Measured Values of Film Thickness 

Measured data is contained in Table III. Pressure readings were 
plotted and the straight line slopes taken to give dp/dz. Fron Equa- 
tion (71), by neglecting Py compared to po, Ce, ean readily be canrputed. 
Similarly .— ean be found from ite definition and the measured values 


of x. Calculations are also inelwied in Table 1:1. 


he 


Lowe att? 
tear 


a Ae 


mar 
Sakari de® 


camry eee Eee 


o aaa rt ned le eh 
rik ce . & oe eo eo © Poe FF Oe’ ae” 
‘e 


~~ + —< —|- ee ed t= oe — 
. ~ | —- oe — a ee oe 























= ; 
os 
o- ‘_- e —— wi « 
» s 
‘s+ -@ ’ al Dekel teddies é 


> = =— a}. cee a of Oa 
*- om i >+S «= =e _—!, 2 





-—-> ° - = 7 : —— 
«Gag? 

_ x, 
_—>- © eae - ‘et » Gia ot 





TABLE Iit 


LIQUID N= 3550 


Water Orifice: 0.10 in. He. 
CPM -093 

Weter Temp: 73.5° F. 

Water Visecsity: 1.953 x 10° ib-see/ft 
Rie 538 

(Yr) 568 1/in. 

Probe well Reading:  .5265 
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445 66 We 52 UW" 1781 5.480 37.3 .5823 
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Liquid Be = 455 


Water Orifice: 
GPM 140 
Water Temp: 10°F. 

Water Viscosity: 2.034 X 107” Ib-sec/rt* 


0.20 in. Hg. 


Me = OP 
( Ay ) i DOL 1/in 
Probe Wall Reading:  .526 
in. 1,0 
inches © belo ate. —_—_—— inches 
a a ir Le a 
ka ing, 
¢) a) ¢) ¢) ¢) 4) 0 ¢) 5015 24.5 13.5 
10 21 oll 27Ti 3 LAG 1. 3.45 524 12.0 624 
12 «14 80 -72 133 1.79 40 4.00 65170 9.0 4.9K 
32 54 63 157 2.13 288  .55 6.52 «5985 7.5 hah 
28 9 1.55 2.77 3.00 4.12 80 9.19 .5195 6.5 3.56 
6L 1.09 2.60 4.22 &.6€2 5.64 1.08 12.35 .520 6.0 3.31 
14h 246 (Ftd) 5.77 862 7.76 HO 16.09 =. 521 5.0 2.76 
2.08 5.07 5.4: 7.22 8.53 Yyobh 1,568 19.03 . 502k 8.6 2.54 
2.85 4.12 7.06 9.20 10.84 11.95 1.95 22.4 5018 4&.e oe. 
5-65 5.05 &.40 10.55 12.43 13.79 2.45 28.2 .9218 be 2. 
5.65 5.22 68.74 11.32 15.20 14.79 2.68 30.7 .5219 4.1 2.26 
4.0 5.85 9.5 W.1 14.20 15.9 2.82 %2.5 .522 4. = 2.20 
11656 15.5 25.4 2.6 34.6 38.6 6.47 74.4 19225 3.5 21.93 
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Liquid Se = 561 


Water Orifice: 0.30 i+. lig. 

CPM LTT 

Water Tenips 1° F. | 

Water Viscosity: 2.010 X 107 

Mee = 51 
( 3/3)” 556 1/in. 

Probe Wall Reading: 5265 

inches ¥,0 below atm. 

Py Po Ps Py Ps 
0 0 O 0 0 
8) .08 2S eli 257 
0 08 050 .28 AL 
@) .08 015 35 (1.0 
0 08 Me) 1.45 2.0 
0 30 lel 1.55 2-2 
9) 65 1.5 2.5 3.35 
15 1.0 2.2 5.5 %:3 
50 1.5 2 95 4.55 5 45 
035 2.25 4.0 59 TR 

1.5 3.0 4.0 6.4 G0 
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28.0 5218 4 2 
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3.4 5222 3.8 
3.9 .5222@ 3% 
45 .6 4825 0 eT 
4 SSe7 83.5 
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Liquid N.. = 747 
water Orifice: 
GPA 0255 
Water Temp: 65° e. 

ai 
Water Viscosity 2.060 X LO 7 Ubesee / et" 


0.50 in. lg. 


= THF 
(2/-*) . 54601 /4n. 
Probe Wall Reading: .5262 
in. B50 
inehes H,© below atm. — — inches 
% 

Py Po Pz, iP Pe p, ap/az ie + ee x. 
© ¢) 0 0 0 0 9) 6) 5015 24.7 13.5 
) 0 03 27 21 8 05 57 $071 19.1 10.43 
Oo 9) 05 52 34 08 1002 «651i 15.1 6.25 
0 6) .03 1% 21.6 1 866152 «64K (CORSE Cee O61 
9) 0 22 72 1.09 1.3 2 $45 527 Qh eee 
0 08 69 1.17 1.98 2.31 .492 5.61 .528 8.2 4.48 
9) -Y 1.10 1.96 2.51 43.31 .725 8.35 .519 7-2 3.93 
0 ce 15 OH 4.05 4.9 21.05 11.95 .5198 6.4 4.50 
“5 ee > ee) a a a ye ne an er) 
&.3 2.9 5.2 715 9.15 10.45 2.15 @4.5 .5@11 5.1 2.79 
2.35 4.35 7-9 10.4 12.45 14.3 2.775 31.6 5215 4.6 2.52 
4.0 6.75 11.2 14.6 17.5 20.55 3.675 41.8 5220 4.2 B29 
5 .B5 9-4 15.1 19.4 23.45 27.3 4.605 52.6 15222 4.io 2.18 
9.0 14.1 621.95 28.2 34.0 %.1 6.075 69.1 .52238 3.4 1.86 
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Liquid Noe = 933 


Water Orifice: 0.60 in. Hg. 

CPM -260 

Water Temp: 69° F. 

Water Viscosity: 2.060 X 107° ib-sec/#t® 
Ree 855 

(ye)? 546 1/4n. 

Probe Wall Reading: .526 


in. 


inches #0 below ata. _ ~ inches 
Film 
PL Po PsP, Pg S/S . Paling xgl0” x 
0 0 0 fe) 9) 9 0 ¢) 5025 23.5 12.83 
.20 0 35 O 65 @ AG 1.8 615 1682 7a 
3 00 6 Sy 15 16 «275 3.45 1.516 10.0 5.46 
JAS 0 6 1.4 2.1 2.9 500 5.70 .517 9.0 %.91 
50 50 18 8.5 SoS 4&2 800 9-11 519 87-0 5.02 
5 > 815 3.5 44 5.2 1.025 121.69 .520 6.0 3.28 
- 1.15 2.5 4.0 5.2 6.2 10085 13.96 58203 5.7 5.4 
~ 1.55 35.15 4.8 6.2 7.25 1.505 17.9% .5208 5.2 2.84 
65 2.05 3.9 5.8 7.4 8.7 1.825 20.8 .5209 3.1 2.78 
1.1 2.60 %85 6.95 8.95 10.25 2.20 25.1 .5210 $5.0 2.75 
1.65 54 6.2 88 10.85 1245 2.75 31.5 15215 4.5 2.86 
2.85 $25 9:06 We we 1711 30 HO Sere 36 28 
5.15 8.2 15.55 17.5 22.3 24.9 4k.6e5 52.6 .5225 3.5 1.92 
9.45 13.45 21.35 27.6 33.0 38.0 6.575 76.0 .5223 5.2 1.75 
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_ 
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Py 
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27 
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4.7 6.2 

7-8 9.9 

11.7 14.8 
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33 
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1.97 

2.70 
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4.37 

5 «50 
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3.78 
3.79 
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9 
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50.1 
64 1 
81.4 
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0521 
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22 12.1 
15 8.25 
10.5 5.78 
@.5 &.é7 
6 3.30 
5 2.75 
4.5 2.48 
4 2.20 
5.5 1.925 
3 1.650 
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